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Abstract 



A class of d-dimensional reaction-diffusion models interpolating continuously 
^2 ' between the diffusion-coagulation and the diffusion-annihilation models is in- 
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\& • troduced. Exact relations among the observables of different models are es- 
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tablished. For the one-dimensional case, it is shown how correlations in the 
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initial state can lead to non-universal amplitudes for time-dependent particles 
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I. INTRODUCTION 



In the recent years, reaction-diffusion problems have stimulated a large body of work in 
many different directions |T]||. The simplest example is provided by the annihilation process 
A + A — > 0, in which the A particles diffuse and react by pairs on contact. 

Usually such kinetic processes are described in terms of macroscopic rate equations giv- 
ing the time evolution for the local averaged concentrations. One assumes that the reaction 
is completely described in terms of the local average densities, and that the reaction in- 
troduces no correlations between the reacting species. This is reminiscent of a mean-field 
like approximation in statistical physics. However, an important aspect of the problem is 
neglected, namely the microscopic fluctuations and it is well known that these fluctuations 
play an important role in low dimensional systems 

Indeed, for the annihilation process A + A — > 0, the rate equation predicts that, in 
the long time regime, the concentration of A will decrease in time as a ~ At -1 while 
a calculation taking into account the local microscopic fluctuations in the particle density 
gives a ~ At a , with a = Min(l, |), where d is the dimensionality of the system. Moreover, 
for a Poissonian initial state, the amplitude A is a universal quantity, which is in particular 
independent on the initial density of the particles. This nonclassical power law behavior is 
called anomalous kinetics. 

Thus to describe correctly these systems, it is crucial to work within a formalism which is 
able to keep track of the fluctuations. This is a difficult task and this is the reason why exact 
analytical results are scarce, but for one- dimensional systems. Complex reaction-diffusion 
systems are present in nature J7|. However, the study of simple models, like the diffusion- 
coagulation process (DC) A + A — > A, or the diffusion- annihilation one (DA) ^4 + ^4^0 
provide a very useful testing ground for new theoretical approaches. 

Different methods have been explicitly developed for the one- dimensional systems (see 
for example |I|j8|-[T3|) . In most of them, one maps the reaction diffusion problem onto a 
different model (quantum chain jl], kinetic Ising model ||, invasion process |I0], etc), which 
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is exactly solvable. A different approach consists in finding the "good observables" for which 
an exact equation of motion can be derived and solved. For the DC model, Doering et al. fnfl 
have shown that the system was best analyzed in terms of the time dependent probability 
S(x,y,t) that the interval between y and x is empty at time t. The particle density a(x,t) 
is simply given by a(x,t) = — ^£(x,y,t)\ x=y . In the continuous limit, £(x,y,t) obeys to a 
simple diffusion equation which can be solved analytically. However, Doering et al. were not 
able to find a similar quantity for the DA model. 

In arbitrary dimension, it was soon recognized, based on dynamical renormalization 



group arguments Jl5|,|T6[, that the DA and DC models belong to the same universality 
class. The power law exponents describing the decay of the particle densities are the same. 
However, for the amplitude A, the situation turns out to be more subtle. Whereas for the 
DC model this amplitude seems to be strictly independent of the initial state, this is no 



longer true for the one-dimensional DA model. Indeed, as noticed by Family and Amar ||T7 
short-ranged initial correlations can change the asymptotic amplitude in a continuous way. 
In their analysis, the DA model is mapped onto a kinetic Ising model at zero temperature || . 
Explicit asymptotic results for the average domain size, average magnetization squared and 
pair-correlation function are derived for arbitrary initial conditions. For the case of an initial 
magnetization ttlq = 0, the results for the DA model with Poissonian initial conditions are 
reproduced. However, for m 7^ 0, the particle density has a nontrivial dependence on the 
initial magnetization. A nonzero value of m means for the DA model that there are some 
correlations among the particles in the initial state. What is the mechanism responsible for 
this dependence and could a similar situation occur also in the DC process for which there 
is no mapping onto a kinetic Ising model? 

To answer this question we shall, in this paper, revisit the problem of universality in 
the DC and DA model. The strategy is to study the field theory which can be associated 
to a reaction-diffusion process, following the method introduced by Doi |L8| and revisited 
by Grassberger et al. |0| and Peliti pt| . We shall consider a class of diffusion-coagulation- 
annihilation models (called a-models), allowing us to go continuously from the DC process 



to the DA one. 

The paper is organized as follows. In Section 2, the way to associate a field theory to 
a reaction-diffusion process is briefly reviewed and the a-models are introduced. Within a 
functional integral formalism, different (i-dimensional a-models are related one to the others. 
In Section 3, the one- dimensional DC model is revisited. It is shown how the equation for 
the time dependent probability £(x,y,t) that an interval between y and x is empty at 
time t can be obtained within the field theoretical framework. In Section 4, we introduce 
the quantity T^ a \x,y,t) which is the natural generalization for the a-model of £{x,y,t). 
Its equation of motion can be derived and solved exactly, leading to an exact expression 
for the asymptotic behavior of the particle density. It is shown explicitly how and when 
the amplitude universality can be violated. In particular, it is proved that for an initial 
state in which particles are pairwise correlated, the violation of amplitude universality is 
only possible for the DA model. Finally, some extensions of this work are discussed in the 
conclusion. 



II. FIELD THEORETICAL APPROACH 

For the sake of completeness and to set up the notations, we briefly sketch the main steps 
necessary to derive the field theory associated to a d-dimensional reaction-diffusion process. 



We follow the method introduced by Doi [18| and revisited by Grassberger et al. [19] and 



Peliti |g| . 



We start from a description of the system in terms of the probabilities P n (x\, . . . ,x n ;t) 
that a configuration with n particles at points x\ , . . . , x n is realized at time t. As we are not 
interested in specifying which particle is where, we choose P n (x\, . . . , x n ; t) to be a symmetric 
function in the variables Xi, . . . , x n . Note that, as the particle number is not conserved, a 
state of the system must be specified by the entire set $(£) = {P n (^)}n=o,...,oo with the 
normalization: 

V / dx l ...dx n P n {xi, ...,x n ;t) = l . 
„ n\ J 



The dynamics is denned by the master equation: 

d t P n (x x ,...,x n ;t) = [HP) n (x u ...,x n ;t) , (1) 



where Ti is an operator acting on the set of the probabilities • 

Let us consider as an example the following DA process. The A particles diffuse and 
annihilate each other with a reaction rate V(\x — y\) depending on the distance. The 
corresponding master equation reads: 

n Q2 

d t P n (xi,...,x n ;t) = D^—^P n (x x , . . . ,x n ;t) (2) 

i=i ax i 

1 If 

-- 22V(\xi - x j \)P n (x x , ...,x n ;t) + -J dydz V{\y - z\)P n+2 (y, z, x x , ...,x n ;t) , 

where D is the diffusion constant. 

We introduce a Fock space representation for the states of the system. Space-valued 
annihilation ip(x) and creation n(x) = ^(x) operators are introduced. The vacuum state 
1 0) is defined by 

ip(x)\0) = 

and the annihilation and creation operators obey bosonic commutation relations: 

ty(x),ir(y)] = 6(x-y) . 

The state specified by the set $ = {P n } n =o,...,oo reads: 

00 1 f 

|$) = 22 — / dx x .. .dx n P n {x x , . . . ,x n )7v(x x ) . . .7r(x n )|0) (3) 

n=0 n - J 

and the statistical average (A) over the state $ of an observable {A n (x x , . . . , x n )} n= o,...,oo 
defined by 

(A) = 2, r / dx x . . . dx n A n (x x , . . . , x n )P n (x x , . . . , x n ) 

n ™ J 

takes the form of a scalar product: 
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where A is the corresponding operator of the Fock space, such that (0\tp(xi) . . . ip(x n )A\&) = 
A n (xi, . . . , x n ) and (| is a projection state given by 

(| = (0|e/ d2 , 

with the property (|7r(x) = (|. In particular, the particle density is given by c(x) = (ifj(x)) 
and the n-point correlation functions by p(xi, . . . , x n ) = (ip(xi) . . . ip(x n )). The normaliza- 
tion of the state is expressed by the condition (|$) = 1. A correlated state, characterized 
by the p-points cumulants: 

cum 9(p) ( x l !■■■■> x p) 



has the form: 



|$) = exp 



00 1 f 

-7 dx 1 ...dx p g(p)(x!, . . . ,x p )(7r(xi) - 1) . . . (n(x p ) - 1) 
p=i P- J 



|0> 



(4) 



It can be shown that the master equation can be written in the Schrodinger-like form: 

d 



dt 



Mt)) = HMt)) 



(5) 



For our annihilation model, the non-Hermitian evolution operator H is given by 



H = Hq + H a 



with 



and 



H = Djdz ir(z)V 2 ip(z) 



H a = - J dzdz' V a (\z- z'\)[k{z)>k{z') - l\^{z)^{z') 



(6) 



(7) 



A second example is provided by the diffusion-coagulation process. The interaction part of 
the evolution operator becomes: 



j dzdz' V c (\z - z'\)[ir(z)ir(z') - n(z)]^(z)^) 



(8) 



If now both coagulation and annihilation reactions are simultaneously allowed, one has: 

H a+C = H a + H c = -J dzdz' [(K + V c )n(z)n(z') - V c tt(z) - V a ]iP(z)^(z') . (9) 

Although not Hermitian, these evolution operators have the property to preserve the state 
normalization: (\H — 0. The evolution equation of an observable A is: 

1{A) = (AH) = ([A,H\) . (10) 



It can be shown ||20|| , that the correlation functions may be expressed by functional 
integrals as: 



. . .ip(x n ,t n )) s ,F = J VipVij) ifj(xi,ti) . . . ip(x n ,t n )e s F [ij(x,0)] . (11) 
The action S appearing in the measure e~ s is related to H as follows: 

S[jj), ij)] = I dtdx ijtp- I dt H[tt = $ + l,ip] . 



F is a functional of the auxiliary field ip at t=0, and is related to the cumulants of the 
initial state by 



F [V>] = exp 



00 1 r - 

— J dx 1 ...dx p </(p)(xi, ...,x p ; 0)^(xi) . ..if>{x p ) 



(12) 



[p=i P' 

Let us now assume that V c is proportional to V a . This leads us to introduce a class of 
model depending on a continuous real parameter a (called the a-models) and defined by the 
following action: 

S a = J dtdz DV 2 i)) + J dtdzdz V(\z - z'\)[i){z)i){z') + a-ip(z)]ip(z)ip(z') . (13) 

The simple coagulation and annihilation reactions correspond respectively to a = 1 and 
a = 2 and will be labeled by the indices c and a in the following. Thus the a-models 
(a G [1,2]) interpolate between the coagulation process and the annihilation one. 

We are now in the position to relate dynamics with different values of a in a very simple 
way. Taking the coagulation reference, eq. (|l3|) gives: 

a 



and, by simply rescaling the functions in the functional integral, we get the following relations 
between correlation functions: 

(ipix!, tx) . . . i[>(x n , t n )) a Fo = a~ n {tp(xi, ti) . . . i[>(x n , t n )) a/a) , (14) 

where: 

F (1/Q) M=F M . 

In particular this means that, if at t = we consider two states |$ a (0)) and |$ c (0)) such 
that: 

(I^On) . . . ^(s„)|$«(o)) = o- n (|^(xi) . . . V(x n )|$ c (o)) , 

for all n, then this relation in conserved for all times, if |$ a ) evolves according to the dynamic 
of the a-model and |<& c ) according to the dynamic of the coagulation one. Moreover, if a set 
of equations among correlation functions, or more generally among operators, is satisfied for 
the coagulation case, then a similar set of relations exists for the a-model, providing that 
one rescales the fields according to ip — > atp. Note that similar relations have been recently 



derived in a different context by Henkel et al. [21 



For the particular case a = 2 and for homogeneous Poissonian initial conditions for both 
|$ c ) and |$ ), we have the following relation between the concentrations: 

c c {t) = 2c a {t) Vt>0 , (15) 

if initially c c (0) = 2c a (0). This result has been already derived by several authors (see for 
example |TT| , p!2|j2T[| ) . Moreover, the relation among the actions S c and S a has been already 



used in previous works [pl||T6|] for Poissonian initial conditions. However, our results are 
valid for all the ct-models and take care explicitly of arbitrary initial condition. In addition, 
from eq. (|HD , we clearly see that quantities involving correlation functions of the same 
order are simply related. But more complicated quantities as for example the interparticle 
distribution function (which involves many correlation functions of different order), are no 
longer simply related. This explain why the asymptotic interparticle distribution function 
between DC and DA model are qualitatively different p2 . 
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III. THE ONE-DIMENSIONAL COAGULATION MODEL REVISITED 



Let us define £(x,y,t), the time dependent probability that the interval (y, x) is empty 
at time t. Our goal is to rederive in the framework of the field theory the equation of motion 
for £(x,y,t) obtained by Doering et al. fl^] . This equation reads: 



d£(x,y,t) (d* a 2 \ 

D \ — + — ) £(x,y,t) . 



(16) 



dt \dx 2 dy 2 

First, we have to identify the Fock space operator corresponding to the probability £(x, y, t). 
This means that we are looking for a time-independent operator E(x,y) such that, given a 
state | the probability that the interval between y and x is empty at time t is given by 



£(x,y,t) = {\E(x,y)mt)) ■ 



(17) 



This operator is: 



E(x,y) 



(18) 



The easiest way to show this consists in using the explicit form of the projection state 
(| = (0|e^-°° dz ^ z ^ in eq. (|T7|). Indeed, for a general state |$), we find: 

£( x , y ) = (0|e ( r + O^)|$) . 



If we consider the state |$') = tt(zi) . . . 7r(^ n )|$) obtained by adding particles in z 1( # 2 , • • • } 2 n 
to the state |$) we obtain: 

f (E(x,y))$ if z n £ (y, x) Vn, 
{E(x,y))& = i 

{ otherwise. 
This defines the probability £(x,y) introduced above. Moreover, we have: 



(19) 



ip(x) 



^E(x,y) 



(20) 



x=y 



in agreement with Doering et al. [14 



The basic commutation relation is: 
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[E{x, y), tt(z)} = -f(z, x, y)E(x, y) , (21) 
where the function / is given by: 

f 1 if y < z < x, 
f{z,x,y) = l (22) 
1 otherwise. 

This commutation relation can be used to prove eq. ([19]) without any reference to the explicit 
form of the projection state. 

We can now derive the equation of motion for £(x, y, t). From eqs. (|SlJT0D and the property 



([POD, we have: 

^(E{x,y)} = (E(x,y){H + H c )) (23) 



= -D(E(x,y)[iP'(x) -iP'(y)]) + / dzdz'V(\z - z'\)(E(x,y)^(z)^(z')) , 
where the domain of integration is 

s = i( z i z ') I z £ (y,v), z ' e (y,x)} 
To proceed, we have to specify the interaction V(£). We choose the following form: 

{V if < £ < (T, 
(24) 
o ife>^ 

and we shall eventually take the a — > limit. For this purpose let us define: 

va = A and va 2 = x ■ 

Taking into account the form of the domain of integration E, the expansion in a of eq 
gives: 

d -{E(x,y)) = {E(x,y)(H + H c )) 

1 



(H (E(x,y)) = (E(x,y)(H + H c )) (25) 



-D(E(x, y)W{x) - ^(y)]> + ^ X (E(x)[iP 2 (x) + ^ 2 (0)]> + 0(a 3 ) 
From the definition of E(x,y) (eq. (|T8|)), we obtain: 

{ & + §i)( E ^v)) = -{E(x,yM'(x)-^'(y)}) + (E(x, y)[^(x) + ^{y)}) . (26) 
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Thus in the limit 

a — ► 0, v — > oo such that i>cr 2 = x — * 2-D , (27) 
we recover the desired equation: 

|(^^)) = ^(^ + ^)(^?/)) + ^ 3 ) ■ (28) 

The particular choice x ~^ 2D has its physical motivation in the relation between a model 
on a discrete lattice and its continuous version p3| . Moreover, our closed equation is valid 



only in the limit Q27D , which corresponds to an infinite pointlike coupling: 

V(\z-z'\)^2\8(z-z'), with A = va -»• oo . (29) 

This fact is not surprising, considering that such a coupling corresponds to an instanta- 
neous coagulation on contact: this is exactly the condition which was used in all the exact 
solutions j8"Hl3 . 



This aspect can also be nicely understood in the renormalization group framework (see 
for example ||). The limit A — > oo corresponds to approaching the non-trivial fixed point 
of the theory and thus to be in the asymptotic long time regime. In consequence, we can 
interpret the result (p8|), valid in the limits (|27j) and (|29"[), as an equation for the asymptotic 
long time regime for a diffusion-coagulation problem, with arbitrary reaction rate. 



IV. GENERALIZATION TO a- REACTIONS 

Following the results obtained in section 2, on the correspondence between the coag- 
ulation model and an a-model, we can now find the quantity corresponding to E for the 
a-model and its equation of motion. 

Let us define the operator: 

T( a \x,y) = e- a fy dzip{z) . (30) 

Thus, T^^x.y) = (T^(x,y)) satisfies the closed equation: 
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at 



T^\x,y,t) = D (^ + £^T^(x,y,t) + 0(a 3 ) . (31) 



To see the physical meaning of this operator, we consider a state |$) = tt(zi) . . . ir(z n )\0). It 
is then easy to show that: 

, , JL f 1 if * £ (y,x), 

(T^\x,y)) 9 = Hai with a t = (32) 
i=i [ 1 — a if Zi G (y, x). 

T^ 2 ^(x, y), which is relevant for the annihilation case, is related to the probability of having 

an even number of particles in the interval between y and x. For a generic vector of the 

form |$) = F[7r]|0), we also have: 

Tto(x,y) = F[l-fV>{z)) , (33) 

where the function is the generalization of fl22"|), i.e., 

fa if y < z < x, 
/<")(*) = \ 

< otherwise. 

Eq. (j3~ID has been solved by Doering e£ a/. [24], with appropriate boundary conditions at 



x = y and x = oo and the initial condition T^ a \x,y,0) = T$ (x,y). The solution of this 
equation is particularly simple if the system is homogeneous. From translational invariance 
T^ a \x,y,t) = T^iyX — y,t) and eq. (]5Tj) becomes: 



^\x,t) = 2D-^T^(x,t) , (34) 



with the conditions: 



(0, t) = 1 , |T< a > (oo, £) | < 1 and T (q) (x, 0) = T (a) I 



x 



The second condition follows from the fact that T^ a \x,y) is a stochastic variable which 
takes its value in the interval [—1, +1]. To take care of the first two boundary conditions, it 
is judicious to define: 

t (-)( x ,t) = -^T^(x,t) , (35) 

whose equation of motion is: 
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f£ (M(M) = 2fl |L iMM) , 

t(«)(0,t) = t^(oo,t) = . 



The solution reads: 



(36) 



t^ a \x,t) = I dz G(x,z,t)t^'(z) 



(«), 



(37) 



where G(x, z, t) is the Green function: 



G(x,z,t) = (SirDt) 



'1/2 



cxp 



(x — z)' 



exp 



(x + z) 



(38) 



8Dt j v 8Dt 

We can now study the behavior of the particle density. From eq. (0), it follows that the 
(homogeneous) density is: 



Ca(t) = (ip{0)) = -a 



-i 



(39) 



x=0 



and then from (j3~5|): 



Ca (t) = a~ l dzt {a \z,t) . 



(40) 



Using eq. fl3T|) we obtain: 



c a (t) = a 1 J dz erf ^ 



(41) 



a- 1 (27r J Dt)- 1 / 2 / d£ e-« 2 e [l - % (a) (2^ 2 ) 
Jo 1 



whose asymptotic behavior for large t is: 



c a (t) ~ (2a)- 1 (2yr J Dt)- 1/2 



i-r w (oo; 



(42) 



This result shows that, besides the trivial dependence on the factor (2a) _1 which follows 
from the rescaling i)j — > a^>, the asymptotic density amplitude depends also on the initial 
value 7q (oo). The study of this quantity for a general initial state is a very difficult task. 
From (|33| ) and (f|) we have: 



T^"- 1 (x) = exp 



00 '-aV> 



(43) 
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where 

PX PX 

9( P )( X ) = J Q dx l ... dx p g {p) (x u ...,x p ) 

and it is difficult to predict the value 7^(oo), without expliciting the g^s. 

Let us then consider a translational invariant state having the following properties: 

g(i)(x)=c, g(2)(xi,x 2 ) = g{\x% - x 2 \) and g(p) = , for p > 3 . (44) 

Such a state can be easily built, as we shall see later. Thus, we have: 

T {a) (x) = exp J -acx + %-g<$ (x) ) , (45) 



with 



If we define 



PX PX 

3(2) (a;) = 2x dz g(z) - 2 dz z g(z) 



PCO POO 

A = 2 dz g(z) and B — l dz z g(z) , 
Jo Jo 

which are finite for short range correlations, one sees that in the limit x — > oo three cases 
have to be distinguished: 



1) 


ac — 


a 2 

Y A>0 


:7< a h 


;oo) 


= 


2) 


ac — 


^A<0 
2 




(oo) 


= oo 


3) 


ac — 




:7< a h 


(oo) 





(46) 



The case 2) is unphysical, because the particle density must be positive. Thus the value of A 
is such that a A < 2c for all a G [1,2], which implies A < c. This constrain is a consequence 
of the assumption g^ = for p > 3 and B < oo. This also means that the case 3) is only 
realizable for a = 2, A = c. From (^) it follows that for an initial state of the type (f44"D, we 
have the same asymptotic density amplitude for all the a-models with 1 < a < 2 and for 
the pure-annihilation case with A < c. But in the pure-annihilation case with A = c, the 
amplitude depends on the value of B. 
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We now show how to built a state of the form ([H]) and understand its physical content. 
Let us take a Poissonian state with density c: 

\P) = e £ -f dz frto-^lo) 



and the state \$p) obtained by the replacement: 

1 



tt(z) — » 7l"(,2) 



< /9 < 1 



(47) 



The normalization = 1 implies that: 



A = 2 dy g(y) = 2(1 - 0)c . 
Jo 

The modification ( f4"?D means that we go from the state \P), where particles are randomly 
distributed, to the state |<&g), where both single particles and pairs are randomly distributed. 
The probability density that a pair have an extension a is ^g(a). If f3 = there are only 
pairs in the state |$ ) • A straightforward calculation gives: 

1^) = exp |(2 - (3)c J dz (tt(z) - 1) + i | rfzrfy ^(z - y)(7r(z) - l)(7r(y) - 1) 

Thus we have a state of the type (0) with density c = (2 — /5)c. For (5 = 0, we have A = c. 
The mean extension of a pair is: 

(2-/9)5 



1 r°° 

° = 7p / da (7#((7) 

2c jo 



2c 



Our a-models, with initial correlations can now be parametrized in terms of (a, (3, a) G 
[1,2] x [0, 1] x [0, oo). The asymptotic density amplitudes are universal for all the values of 
the parameters, except on the line (a = 2, j3 = 0, a) where the amplitude is proportional 



to (1 



-4C(T X 



. In this case, the asymptotic density amplitude tends to zero if cr — ^ 0. The 



physical reason is that when a — 0, the two particles of a pair are at the same place and 
immediately annihilate. This situation cannot occur if coagulation is present, because single 
isolated particles will remain. 

It would be interesting to see how a more general initial state would affect the amplitudes. 



Note finally, that the relations fl46|) , characterizing the initial state, are preserved by the 
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dynamical evolution. Indeed, every evolved state can be considered as a new initial state 
evolving towards the same asymptotic solution. 

V. CONCLUSIONS 

We have studied a class of d- dimensional diffusion-reaction models interpolating contin- 
uously between the diffusion-coagulation and the diffusion-annihilation models. The field 
theoretical approach used leads to exact relations between the observables of these different 
models. In one dimension, it was shown how correlations in the initial state can lead to a 
violation of the universality of the amplitude for the DC models. 

Several extension of the present work concerning one-dimensional models with reversible 
diffusion-reaction systems, or the presence of fronts in inhomogeneous systems are under 
investigation. Moreover, one may expect to construct "good observables" in a?-dimensional 
systems for which a closed equation of motion can be derived. 
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